Logarithms of a in QED bound states from the renormalization group 

Aneesh V. Manohar* and Iain W. Stewart^ 
Department of Physics, University of California at San Diego, 
9500 Gilman Drive, La Jolla, CA 92099 

The velocity renormalization group is used to determine In a contributions to QED bound state 
energies. The leading order anomalous dimension for the potential gives the a In a Bethe logarithm 
in the Lamb shift. The next-to-leading order anomalous dimension determines the a 6 In a, aJ ln 2 a, 
and a 8 In 3 a corrections to the energy. These are used to obtain the a 8 ln 3 a Lamb shift and a 7 In 2 a 
hyperfine splitting for Hydrogen, muonium and positronium, as well as the a 2 In a and a 3 ln 2 a 
corrections to the ortho- and para-positronium lifetimes. 



The energies of QED bound states such as positro- 
nium, muonium or Hydrogen depend on In a For 
example, the Lamb shift g{ of Hydrogen contains the 
famous Bethe logarithm of order a 5 In a. In relativistic 
scattering at high energy E 3> m, logarithms of E can be 
determined from renormalization group equations. The 
leading logarithmic series of the form (a\nm/ E) n can be 
summed by integrating the one-loop anomalous dimen- 
sion, the subleading series a(ahim/ E) n by integrating 
the two-loop anomalous dimension, and so on. In this 
letter, we show how to predict logarithms of a for non- 
relativistic QED bound states using the velocity renor- 
malization group (VRG) ||. Our approach can also be 
used for other non-relativistic systems. 

The leading order (LO) anomalous dimension gener- 
ates the Bethe logarithm of order a 5 In a for the Lamb 
shift, and is the only term in the series. Integrating the 
next-to-leading order (NLO) anomalous dimension deter- 
mines the a 6 In a, a 7 In 2 a and a 8 In 3 a terms, and simul- 
taneously gives the Lamb shift and hyperhne splitting 
for Hydrogen, muonium and positronium, as well as the 
decay widths for ortho- and para-positronium. Here we 
derive all these terms except the a 7 hi 2 a Lamb shift and 
a 6 In a Lamb shift and hyperhne splitting, which depend 
on terms in the NLO anomalous dimension for which 
values are not presented here. The a 8 In 3 a term for Hy- 
drogen has been the subject of a recent debate in the 
literature. An analytic calculation by Karshenboim ||] 
and a numerical calculation by Goidenko et al. || agree 
with each other, but disagree with numerical calculations 
by Mallampalli and Sapirstein || and by Yerokhin |t]] 
(which agree). The renormalization group answer agrees 
with the analytic a 8 ln 3 a result of Karshenboim, and we 
will comment on the disagreement with Refs. [g|Jj. The 
a 8 ln 3 a positronium Lamb shift is a new result. 

Our calculation makes use of the non-relativistic ef- 
fective theory for QED @ (NRQED), formulated as in 
Refs. (Hi- We will consider the interaction of two parti- 
cles of mass mi and 777,2, and charge — e and Ze respec- 
tively. The effective theory has a subtraction velocity v, 
rather than the usual subtraction scale /x of dimensional 



rcgularization. QED is matched onto the effective the- 
ory at v — 1, and the effective Lagrangian is scaled down 
to v = Za, the typical velocity in the bound state, us- 
ing the VRG. (Terms oc ln(m2/mi) are suppressed for 
7722 3* 777i and are neglected.) Finally, the energy is 
computed from matrix elements of the Lagrangian. In 
this approach, all logarithms arise from renormalization 
group running, since matrix elements of operators renor- 
malized at v = Za and matching coefficients at v = 1 
contain no large logarithms. This should be compared 
with the traditional approach where In a terms are de- 
termined by examining integrals for the matrix elements. 
Determining the logarithms using the VRG gives them a 
universal and simple description. 

The interaction potential has an expansion in powers 
of the velocity v, 

V = V<-V + W°> + V™ + V<$ + V® + . . . , (i) 

where V^ n ' ~ v n . The lowest order interaction is the 
Coulomb potential, V^ -1 ), of order a/ v. [In momentum 
space V^ 1 ^ ~ 1/k 2 ~ v~ 2 . Shifting the power by one 
makes the power counting simpler.] The potentials in 
the center of mass frame for scattering of particles with 
momenta ±p to ±p', with k = p' p are 
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V^ = (U 3 + U 3s S 2 ) |k| + 

where the total spin S = (cri + <t 2 )/2. could have a 
term of the form p 2 p' /k 2 , but it is more convenient to 
rewrite this as (p 4 + p' )/(2k 2 ) plus terms that vanish 
on-shell. In Eq. (^J) all momentum dependence is explicit, 
and only the vector coefficient Ua depends on spin. The 
odd terms, V^ 2n+1 \ have coefficients of order a from tree 
level matching, while the even terms, V^ 2k \ are first gen- 
erated at one-loop. Matching on-shell at v = 1: 
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The reduced mass 



mim2J m\m, 2 
= rti\m<2,l(m,\ + m 2 ). Terms 
not needed have been omitted. Let U 2 + s = U2 + 
U S S 2 . Positronium has additional contributions from 
one-photon annihilation, U 2+S (l) = iraS 2 /m 2 , and from 
two-photon and three-photon annihilation graphs which 
give the imaginary terms (where m e is the electron mass): 

_ W(S 2 -2) iroV^ 
U2+s ~ m 2 e ' U2+s ~ 9irm 2 e ' 1 j 

Evaluating the matrix element of a potential of order 
a r v s gives a contribution to the energy at order a r+s+2 . 
The iteration of two potentials of order a Tl v Sl and a T2 v S2 
gives a potential of order a Tl+T2 v Sl+S2 . Thus, V^ 1 ', 
y~(°) ; y( 2 ), and contribute starting at order a 2 , 
a 4 , a 4 , a 6 , and a 6 respectively. The products y(°)y(°), 
V {Q) V {1 \ and V {1) V {1) contribute starting at order a 6 . 
In evaluating matrix elements a ~ w, y^ 1 ) - 1, so the 
Coulomb potential cannot be treated as a perturbation. 

The LO anomalous dimension for a potential coeffi- 
cient of order a r will be defined by terms of order a r+1 , 
the NLO anomalous dimension by terms of order a r+2 , 
etc., rather than by the conventional definition in terms 
of the number of loops. In bound states where the typ- 
ical momentum is smaller than the electron mass, the 
Coulomb potential has no anomalous dimension. The 
first terms that have an anomalous dimension are 
and The LO anomalous dimension for V^ ' and 

generate the LO series in the energy a 4 (alna)". 
Their NLO anomalous dimensions generate the NLO se- 
ries a 5 (a In a) n . Matrix elements of or or prod- 
ucts of V(°) and first contribute at order a 6 (a In a)" - , 
the same order as the NNLO anomalous dimensions for 
y(°) and V {1 \ and are not needed for our analysis. How- 
ever, the potential mixes into at NLO, and is 
therefore necessary for solving the VRG equations. 

The anomalous dimensions in the effective theory can 
be generated by soft, potential and ultrasoft loops, with 
energy-momentum of order (Ar ~ mi>,k ~ mv), (k ~ 
mw 2 ,k ~ mv) and (k ~ mw 2 ,k ~ mv 2 ) respectively. 
There are no potential anomalous dimensions at LO. In 
NRQED soft vertices are labelled by a > 1, and a loop 
graph with two soft vertices is of order v ai+a2 ~ 1 ||. The 
one-loop soft anomalous dimension is non-zero only for 
y-( 2 ™-i) with n > 1. At LO only the running of is 
needed (from Fig. |l|a with two a = 1 vertices): 

, ■ (5) 

dv Sm\m 2 

where the other coefficients in V' 1 ) have zero soft anoma- 
lous dimensions. The LO ultrasoft anomalous dimension 






FIG. 1. Graphs contributing to the soft (a) and ultrasoft 
(b,c,...) anomalous dimensions at leading order. The poten- 
tial is denoted by ®, and one sums over all possible ultrasoft 
exchanges (including wavefunction renormalization) . 



is independent of the momentum structure of the poten- 
tial (from Fig. 0b, etc.): 

vfv(p, p') = ^ ( — + — V k 2 V ( P , p') . (6) 

dv Sir \mi m 2 J 

At LO, the VRG equations for v'- ' 1 ^ are: 

dU k 
v—— = 0, 
dv 

dU 2 2a ( 1 Z \ 2 rr UZ 2 a 2 



v-r = it- — + — u c + , 

dv Sir \mi m 2 J Sm\m 2 
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= t" — + — U k + ll U c + l2 U 2 , (7) 
dv Sir \mi m 2 J 

obtained by summing Eqs. (|) and (§). The last equa- 
tion also has additional spin-independent contributions 
denoted by 71,2, which are not needed here. The U 2 in- 
tegration is trivial, 



U 2 {v) = l0 U c \nv + U 2 {l) , 
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and U 2 {1) has no large logarithm. The matrix element of 
the U 2 potential is (mnZa) 3 /irn 3 for the nS state (and 
zero for L ^ states). With the coefficient U 2 (v = Za) 
from Eq. (Bh this gives the logarithmic energy shift 
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which is the well-known Bethe logarithm in the Lamb 
shift, and is valid for Hydrogen, muonium and positro- 
nium. Equation (||) has no imaginary part or spin- 
dependence, so there is no contribution to the decay 
width or hyperfine splitting at this order. Equation (|^) 
has been computed before using an effective field the- 
ory po[ . The VRG method makes it clear that the LO 
series a 4 (a In a)™ has only a single term — the anomalous 
dimension for U 2 depends on a and U c , both of which do 
not run, so integrating the VRG equation produces only 
a hiv term. Below the electron mass, QED is not very 
efficient at generating logarithms. 

At NLO, one needs the anomalous dimension of to 
order a 4 , and of to order a 3 . The possible terms that 
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FIG. 2. Examples of graphs contributing to the NLO 
anomalous dimensions. The Cg> denotes insertions of terms 
in the potential such as U c , U2, etc. 



can contribute are determined by using the ^-counting 
formula in Ref. ||. The NLO anomalous dimension of 
has the form vdUk/dv — jsUc, where 73 is gener- 
ated by Fig. (||)b iterated with U c . The NLO anomalous 
dimension for has several contributions from poten- 
tial loops (see Fig. ||a): U 3 with a v 4 insertion, U 2 V^ 
with a v 2 insertion, U C V^V^\ U%V®>, U C V^ with a 
v 2 insertion, U C V {2 \ and V {1) V {Q) . There is no ultrasoft 
anomalous dimension for V"W a t NLO, but there is a 
soft contribution shown in Fig. |^b with two a = 1 inter- 
actions, which has the form p s Z 3 a 3 /mim 2 . 73 and p s do 
not contribute to our results. The only NLO anomalous 
dimension needed is: 
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where U 2+s = U 2 + U S S 2 and 
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The potential does not mix into when the ba- 
sis is chosen to avoid the p 2 p /2 /k 2 term. Equation (||) 
depends on V^ 2 \ so we also require its LO VRG equa- 
tion in Eq. ((?]) . The ellipses in Eq. (||) denote terms other 
than U3 in the potential (e.g. U3 S ), which contribute 
to the running of U 2 . However, they do not have a LO 
anomalous dimension, and thus contribute with one less 
logarithm and are not required for our analysis. 

Integrating Eq. (||) generates an infinite series of terms, 
since it is non-linear in U 2+s , including a 2 (a\nv) n terms 
with n = 1,2,3. Thus, the NLO a 2 (alnv) n series ter- 
minates after three terms. Integrating the U C U 2 term in 
Eq. (^) gives these three possible terms, depending on 
whether each U 2 is replaced by the joU c \ni/ or U 2 (l) 
term in Eq. m). 
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FIG. 3. Fig. (a) is the a 8 ln 3 « contribution to the Lamb 
shift computed by Karshenboim. Fig. (b) is the additional 
contribution of Yerokhin. The dashed lines are Coulomb ex- 
changes. 



The simplest term is the hi 3 u term in U 2 + S (v): 
\llPc22U 3 c {l) In 3 !/, 
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whose tree-level matrix element determines the a 8 In 3 a 
Lamb shift for the nS state 
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The a 8 In 3 a term is spin-independent and has no imagi- 
nary part, so it does not contribute to the hyperfine split- 
ting or the positronium lifetime. The structure 7gp C 22 in 
Eq. (]ll]) is the same as the matrix element in Fig. ||a 
computed by Karshenboim and our results reduces 
to his for Hydrogen (m 2 — > 00). For positronium, we 
find AE = 3m e a 8 ln 3 a/(87r 2 n 3 ). In the effective theory, 
the matrix element (for the log terms) becomes the prod- 
ucts of the graphs in Fig. ||a and Fig. [|. Yerokhin has 
an additional contribution Fig. ^b, which has the struc- 
ture 7g times the U 2 matrix element. While the loop 
graph with two U 2 's is linearly divergent, it is finite in 
dimensional regularization and does not contribute to the 
anomalous dimension or to the a 8 In 3 a term. 

The piece of the In 2 v term that depends on the match- 
ing value for U 2 + s at v = 1 is 



70Pc22 U 2 (l) \U 2 (1) + U s (l)S 2 ] In 2 v. 
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The remaining terms depend on the matching values for 
U c and Uk, are spin-independent, and only contribute 
to the Lamb shift. The matching value is U s (l) — 
AnZ pip 2 a/ (3mim 2 ) where pi are the magnetic moments 
in units of ej (2m,i). Using Eq. ( |l2| ) and taking the differ- 
ence between the energies for S 2 = 2 and S 2 = gives the 
a 7 In 2 a hyperfine splitting for muonium and Hydrogen 
for the nS state: 
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in agreement with Refs. |4],|11|]. The matching value for 
positronium, U s (l) = 77rcv/(3?77 2 ), and p\ = p 2 = 1 gives 
the hyperfine splitting for the nS state: 
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(13) 
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in agreement with Refs. |M l2Ul3| . 

The positronium width is T = — 2lmE. The match- 
ing coefficients U 2 + s have imaginary parts from the two- 
and three-photon annihilation graphs, of order a 2 and 
a 3 respectively. They are the lowest order contributions 
to the imaginary part, so we can use Eq. for the 
widths, even though real parts of the same order have 
been neglected: 



Solving this equation gives 



■77- = 7oPc22 U C {1) 2 In 2 v = --^-a 3 ln 2 a, 

1 2,71 



(14) 



for the ortho- and para-positronium widths in agreement 
with Ref. [Q. The a 7 In 2 a Lamb shift depends on 7i i2 , 
and will be discussed elsewhere. 

The In v term that depends on the matching value of 
U 2 + s is 
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where the L^'s are evaluated at v = 1. Using Eqs. JIo| ) 
and (0) , the imaginary part gives the positronium width: 
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in agreement with Ref. [Q . The a 6 In a Lamb shift and 
hyperfine splitting have contributions from V^^v = 1), 
73 and p s . 

There are two infinite series of logarithmic terms that 
are easily identified. Neither of them gives the complete 
contribution at a given order, but they do show that there 
are logarithmic terms of arbitrarily high order. The ul- 
trasoft anomalous dimension Eq. (^|) generates the po- 
tential 
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where the exponential is related to the Sudakov form 
factor. In position space, the modified Coulomb potential 
is 
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where Erf is the error function. This gives a series in the 
energy of the form a b In Za(a 3 In Za)™/ 2 , n > 0. The 
second infinite series of logarithmic terms is obtained by 
integrating the VRG equation for U 2 retaining only the 
f7 2 terms in Eqs. ([?]) and (||), 
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which has an expansion of the form a 2 lna(a 3 In 2 a) n . 
The first two terms are the a 5 In a and a 8 In 3 a Lamb 
shifts. The difference between the exact expression and 
the first two terms is less than 1 Hz in the Lamb shift. 

We have derived logarithmic terms in QED bound 
states using the VRG. The results are in agreement with 
results computed previously using other methods. The 
computation of the a 8 In 3 a Hydrogen Lamb shift by the 
VRG supports the value computed by Karshenboim [||. 
We have also computed the a 8 hi 3 a positronium Lamb 
shift. The VRG approach makes the universal nature 
of the In a terms clear, and is also an efficient way of 
computing these terms. The computations in this paper 
support a key feature of the VRG method proposed in 
Ref. S: the simultaneous running of soft and ultrasoft 
effects from m using a subtraction velocity. We have also 
identified two infinite series of logarithms in QED. 
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